Spin-polarized surface states close to adatoms on Cu(lll) 
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We present a theoretical study of surface states close to 3d transition metal adatoms (Cr, Mn, Fe, 
Co, Ni and Cu) on a Cu(lll) surface in terms of an embedding technique using the fully relativistic 
Korringa-Kohn-Rostoker method. For each of the adatoms we found resonances in the s-like states 
to be attributed to a localization of the surface states in the presence of an impurity. We studied 
the change of the s-like densities of states in the vicinity of the surface state band-edge due to 
scattering effects mediated via the adatom's rf-orbitals. The obtained results show that a magnetic 
impurity causes spin-polarization of the surface states. In particular, the long-range oscillations of 
the spin-polarized s-like density of states around an Fe adatom are demonstrated. 

PACS numbers: 73.20.At, 73.20.Hb, 73.22.-f 
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INTRODUCTION 

Following the first experimental observation of a sur- 
face band at Cu(lll) in terms of angle-resolved photoe- 
mission (ARPES) electrons at closed packed surfaces 
of noble metals have been at the center of much exper- 
imental 0, S S S E and theoretical H H E3 at- 
tention. For a pristine surface the energies of Shock- 
ley surface states lie in the 'gap' around the L-point 
of the bulk Brillouin Zone, their wavefunctions being 
confined to the surface. The corresponding dispersion 
relations ave been determined by ARPES and were 
found to be two-dimensional free-electron like parabo- 
las 0, ■ One of the highly interesting features of this 
phenomenon is the response to perturbations causes by 
placing, e.g., adatoms on the surface. As to be expected 
such a response is characterized by long range Friedel 
like charge oscillations governed by the dispersion rela- 
tion of the two-dimensional surface electron gas. The 
existence of long range interactions between adatoms on 
surfaces supporting surface states can lead to the for- 
mation of an atomic superlattice as recently shown both 
experimentally and theoretically [ljj. Until recently 
STM studies of atoms on well defined Cu, Ag and Au 
(111) surfaces imaged only the charge distribution of the 
surface electrons y, 4j |5| . Remarkable developments in 
spin polarized STM ll|, however, are expected to detect 
spatial variations in the magnetic density, which in turn 
might provide an understanding of the magnetic behav- 
ior of this kind of systems [12, EH Q E| • Evidently, this 
also opens up the possibility of designing magnetic nano- 
structures for both scientific and technological purposes. 

In accordance with a theoretical prediction by Borisov 
et al. 0, by using STM at 7 K, Limot et al. \vj& 
very recently reported on an adatom-induced localization 



of the surface states electrons on Cu(lll) and Ag(lll) 
surfaces. The existence of such resonances for a single 
Cu adatom was also shown theoretically by Olsson et al. 
|l8l | in terms of a parameter-free pseudopotential method. 
The appearance of a peak in the density of states (DOS) 
just below the surface state band can be attributed to a 
theorem by Simon claiming the existence of a bound state 
for any attractive potential in two dimensions |l8lll9| . It 
was demonstrated in Ref. |l7j | by comparing results for Co 
and Cu adatoms on Cu(lll) that the type of adsorbate 
influences the shape of the adatom-induced resonance. 
For this very reason we performed a systematic study for 
a series of 3d impurities (Cr, Mn, Fe, Co, Ni and Cu) on 
Cu(lll) in order to identify spin-dependent features of 
the occurring resonance. 



METHOD OF CALCULATIONS 

Within multiple scattering theory the electronic struc- 
ture of an ensemble of non-overlapping potentials is de- 
scribed by the so-called scattering path operator (SPO) 
matrix (for more details see, e.g., Ref. |2fl). The SPO 
matrix tq that refers to a finite cluster C embedded into 
a host system can be obtained from the following Dyson 
equation |2lj . 



t c (E) = r h {E) [I - (t-\E) - tc l {E))T h {E)]- 



(1) 



where th{E) and Th{E) denote the single-site scattering 
matrix and the SPO matrix for the unperturbed host 
confined to sites in C, respectively, while tc denotes the 
single-site scattering matrices of the embedded atoms. 
Note, that Eq. (JJJ takes into account all scattering events 
both inside and outside the cluster. 
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First, a fully self-consistent (SC) calculation is per- 
formed for the Cu(lll) surface by means of the screened 
Korringa-Kohn-Rostoker (SKKR) method Then 
for single adatoms placed on top of Cu(lll) the multiple 
scattering problem is solved self-consistently in terms of 
the embedding method discussed in details in Ref. [2~l| . 
The selfconsistent calculations for both the semi-infinite 
Cu(lll) host and the adatoms were performed using the 
atomic sphere approximation (ASA) and the local spin- 
density approximation in the parameterization of Vosko 
et al. 23| . Due to the fully relativistic treatment applied 
the orientation of the effective magnetic field had to be 
specified: it was chosen to point along the z axis (nor- 
mal to the surface). In order to evaluate the inevitable 
Brillouin-zone integrations, in the selfconsistent calcula- 
tions 70 fc|| -points were used in the irreducible wedge of 
the surface Brillouin-zone (SBZ). For the calculation of 
the i-matrices and for the multipole expansion of the 
charge densities, necessary to evaluate the Madelung po- 
tentials, a cut-off of imax = 2 was assumed. Energy 
integrations were performed by sampling 16 points on a 
semicircular contour in the complex energy plane accord- 
ing to an asymmetric Gaussian quadrature. The densities 
of states were calculated parallel to the real energy axis 
with an imaginary part of 5 mRyd by sampling about 
2200 fc|| points within the irreducible wedge of the SBZ. 
In the present calculations no attempt was made to in- 
clude surface relaxations: the geometry was taken to be 
identical to an ideal Cu bulk fee lattice (lattice constant 
a = 3.614 A). 



RESULTS AND DISCUSSION 

In order to determine the dispersion relation and the 
effective mass of the surface electrons Block- spectral func- 
tions (BSF) ji^l were evaluated in terms of the SKKR 
method for k-points between the V and K points in the 
fcc(lll) SBZ close to the Fermi energy. It should be 
recalled that the dispersion relation of the surface state 
band can be defined by the position of the maxima in 
the BSF. In agreement with experiments the calculated 
dispersion relation is a free-electron like parabola as in- 
dicated in Fig. n The bottom of the calculated sur- 
face state band is only about 0.3 eV below the Fermi 
energy which is a bit smaller than the experimental val- 
ues (0.4 eV) 0, El Also in good agreement with the 
experimental data Q is the corresponding effective mass 
m* = 0.37 m e as obtained using an appropriate fitting 
procedure. In Fig.^also the density of the surface states 
in the vicinity of the band-edge is displayed. Shown is 
the s-like DOS of the first vacuum (empty sphere) layer 
as integrated over a sphere of a radius of 0.15/a.u. cen- 
tered around f of the SBZ. Note that the corresponding- 
DOS at the substrate layers decays exponentially with 
increasing distance from the surface. As our calcula- 



tions do not include structural defects (e.g., steps) at 
the surface, electron-electron interaction beyond the den- 
sity functional theory or electron-phonon interaction, the 
rather broad onset of the surface-state band is a direct 
consequence of using an imaginary part of 5 mRyd for the 
energy when calculating the DOS, see the solid line on 
the right of Fig.^ In order to justify this argument, only 
for this particular case, we also performed calculations 
with imaginary parts of 2, 1 and 0.5 mRyd displayed in 
Fig. in terms of dash-dotted, dashed and dotted lines, 
respectively. As can be seen, by decreasing the imag- 
inary part of the energy the onset of the surface-band 
approaches a step-like behavior. The experimental onset 
(~ 30 meV |Tjj ) is fairly well recovered in case of 1 mRyd 
for the imaginary part. It should be noted that, in or- 
der to smooth out spurious oscillations in the calculated 
density of the surface states, in this case a sampling of 
more then 50 000 ku points within the irreducible wedge 
of the SBZ (about 2000 fc|| points for \k\\\ < 0.15) was 
needed, see also a recent theoretical STM study of Hofer 
and Garcia-Lekue j2^|. Fortunately, however, applying 
a Lorentzian broadening of 5 mRyd (~ 70 meV) for the 
DOS turned out to be sufficient to resolve the adatom 
induced surface states to be studied in this paper. 
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FIG. 1: Left panel: Positions of maxima in the Bloch spec- 
tral functions near to the T point of the SBZ (dots). The solid 
line refers to a parabolic fit. It should be noted that only the 
first third of the SBZ is displayed (K » 0.65). Right panel: 
Density of the surface states < 0.15) close to a clean 

Cu(lll) surface. The different values for the imaginary part 
of the energy used are displayed explicitely. The Fermi energy 
and the bottom of the surface state band are indicated by a 
dashed and a dotted horizontal lines, respectively. 

In Fig. [3 the calculated s-DOS and total DOS (insets) 
for the chosen series of adatoms are presented as projec- 
tions with respect to the two spin directions. As is also 
evident from the spin-split peaks in the total DOS, domi- 
nated mainly by <i— like contributions, the Cr, Mn, Fe, Co 
and Ni adatoms on Cu(lll) were found to be magnetic 
with spin-moments of S§ r = 4.21 /i^, S|j n = 4.39 /is, 
= 3.27 (i B , S^ Q = 2.02 [i b and = 0.51 /is, while 
the Cu adatom is non-magnetic. It should be noted that 
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within a relativistic description the electronic spin is not 
a constant of motion. In the present cases, however, char- 
acterized by large exchange splittings and weak spin-orbit 
interactions it is quite illustrative to view the two spin- 
channels separately by making use of the approximation 
discussed in Ref. l2ol. 
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FIG. 2: Calculated spin-projected s-like density of states at 
the impurity positions. Insets: Spin-projected total density 
of states of the adatoms. In each entry, the majority and the 
minority spin DOS's are depicted by dotted and solid lines, 
respectively. 

Concentrating first on the case of the Cu adatom, quite 
a broad peak is found just below the bottom of the sur- 
face state band (E —0.025 Ryd). The appearance 
of this peak indicates that adatoms act as an effective 
attractive potential for surface state electrons inducing 
thus a certain degree of localization. The widths and po- 
sitions with respect to the surface state band edge of the 
adatom-induced peak agrees well with the one obtained 
in Ref. 0] . The calculated resonance width being also in 
agreement with Ref. [lij is, however, approximately three 
times larger than experimentally measured [TtI ] . This dis- 
agreement with the experiment can most probably be at- 
tributed again to the rather large Lorentzian broadening 
of the calculated DOS. 



The fairly narrow d-band of the Cu adatom lies well 
below the 2D surface band (Ed ~ —0.1 Ryd). It can 
therefore be assumed that the adatom-induced resonance 
at E w —0.025 Ryd is hardly influenced by d-like states. 
This allows us to consider the Cu adatom as a refer- 
ence when identifying effects of <i-states on the adatom- 
induced resonance in the case of other impurities. In line 
with Ref. 0| , a peak in the s-DOS of the Cu adatom 
can be observed just at the peak position of the d-band. 
This is a direct consequence of the hybridization between 
s and d type atomic orbitals as illustrated in Fig. [21 us- 
ing group theoretical arguments. It is important to note 
that in order to make use of a well-defined non- relativistic 
classification of the densities of states in terms of sym- 
metry adapted spherical harmonics, only in this case we 
"switched off" the spin-orbit coupling by applying the 
so-called scaling scheme proposed by Ebert et al. |27| . 
In the upper panel of this figure the rf-DOS is decom- 
posed into three components: two of them correspond- 
ing to the two-dimensional irreducible representation (E) 
of the C$ v point-group, namely, according to the follow- 
ing sets of symmetrized basis functions {d xy , d x 2_ y 2} and 
{d XZl d yz }, and one to a one-dimensional representation 
(Ai), namely, with respect to {d z 2}. Clearly enough, the 
s-states that correspond to A\ symmetry can hybridize 
only with d z 2 -states as is apparent from the line-shapes 
of the corresponding DOS's in Fig. [31 
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FIG. 3: Upper panel: Calculated d-DOS of a single Cu 
impurity on a Cu(lll) surface as partitioned into d xz , d yz , 
d xy , d z 2 and d x 2_ y 2 like contributions. Lower panel: s- 
DOS at the same site. The DOS's presented in this figure 
correspond to one of the spin-projections. 

From Fig. |3 the appearance of the adatom-induced lo- 
calized state is also obvious in the case of the magnetic 
impurities as only a slight variation of its position in en- 
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ergy occurs. This indicates that each kind of impurity 
acts in a similar manner as an attractive potential well 
for the 2D electron gas. As can be seen, however, in 
the case of magnetic adatoms the shape of the adatom- 
induced resonance is different for the two spin projec- 
tions. This observation can be well explained in terms of 
the energetic position of the corresponding ri-bands of the 
impurities. For Mn, Fe and also for Co the majority ri- 
band lies reasonably deep below the adatom-induced res- 
onance. Therefore, the line shape of one (say majority) 
spin-projected s-DOS around E = —0.025 Ryd is prac- 
tically unchanged as compared to the Cu case. Within 
the energy range displayed, an s-d hybridized peak in 
the majority s-DOS is seen for Co around -0.09 Ryd. 
An analysis in terms of symmetry-adapted spherical har- 
monics (basis functions) as discussed in the case of the 
Cu adatom applies also in this case. For Cr and Ni the 
position of the majority d-band is just at the bottom 
of surface band. The adatom-induced resonance appears 
therefore for Ni just as a small shoulder on the upper side 
of the s-d hybridized peak, while for Cr the resonance can 
hardly be traced at all. 

For all magnetic adatoms investigated the minority ri- 
band overlaps with the 2D surface state band. In the 
series Cr to Ni its position moves downwards in energy, 
its width decreases monotonously. For Mn and Fe a mod- 
erately well-developed peak can be seen at the position 
of the very sharp ri-band, while for Co and Ni just a 
small 'hump' in the minority s-DOS is visible. The oc- 
currence of these structures is of similar origin as the 
s-d hybridized peaks in the majority spin channel: the 
continuum-states experience resonant scattering due to 
an overlap with an adatom's ri orbitals giving thus rise 
to well-known virtual bound state resonances. As can 
be shown in terms of a non-interacting Anderson model, 
this effect is proportional to the DOS of the continuum 
band, which decreases if the position of the ri-band moves 
towards the bottom of the surface band. This explains 
qualitatively the trend for the width of the minority ri- 
band. It is, however, apparent that the minority spin 
adatom-induced resonance peak (E w —0.025 Ryd) is 
considerably larger for Co and Ni as for Cu, a fact which 
might be attributed to a strengthening of the localization 
of the 2D surface electrons due to d— like virtual bound 
states very close in energy. 

In Fig.0]the s-DOS for both spin directions is shown 
for a site at selected distances measured with respect to 
an Fe adatom in the same plane (parallel to the surface) . 
It can be observed that the line shapes related to the 
adatom-induced resonance and the minority spin reso- 
nant scattering vanish in fact at a radius of about two in 
units of the 2D (in plane) lattice constant. Beyond this 
radius a lincshape characteristic for the 2D surface-state 
band appears to evolve, however, with a broader onset 
than for the pristine Cu(lll) surface. This broadening 
can again be attributed to the interaction (overlap) be- 



tween the surface states and the ri-states of an adatom. 
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FIG. 4: Calculated minority spin (upper panel) and majority 
spin (lower panel) s-DOS 's at a site in the same plane as a 
single Fe impurity on a Cu(lll) surface with respect to the 
distance between these two sites. The distance (in units of the 
two-dimensional lattice constant) is indicated in the legend 
of the lower panel. 



At a distance of seven 2D lattice constants (~ 17.9 A) 
the s-DOS is practically the same as calculated for the 
clean Cu(lll) surface. A long-ranged oscillatory behav- 
ior of the densities of states can, however, still be re- 
solved. In Fig. El the DOS at a selected energy, namely, 
34 mRyd above the bottom of the surface-state band 
is displayed for both spin channels as a function of the 
distance from an Fe adatom. A simple estimate of the 
wavelength of the 2D Friedel oscillation, 



hV2m*E 



(2) 



gives A ~ 15 A k 6a2D that can be read off from 
Fig. [S] As shown in the inset of Fig. [31 a magnetic im- 
purity also induces long-range oscillations in the magne- 
tization density of states (MDOS, defined as the differ- 
ence of the spin-projected DOS's), with the same period. 
Clearly, these oscillations in the MDOS lead to a long- 
range RKKY interaction on noble metal (111) surfaces 
as discussed, e.g., in Ref. [Io| . 
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adatoms (Cr, Mn, Fe, Co, Ni and Cu) on a Cu(lll) sur- 
face. In agreement with previous theoretical 
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FIG. 5: Calculated minority spin (solid line) and majority 
spin (dotted line) s-DOS's at E - E F = 0.012 Ryd of a site 
positioned in the same plane but at a certain distance from a 
single Fe impurity on a Cu(lll) surface. Inset: difference of 
the spin projected s-DOS's. 



CONCLUSION 

We presented a series of calculations for spin- 
dependent surface states close to 3d transition metal 

[lj and 

experimental |17| studies the appearance of an adatom- 
induced resonance just below the surface state band edge 
is shown. The occurrence of an additional peak in the s- 
DOS caused by the adatom's d z 2 -states was found be- 
low the adatom-induced resonance not only for a Cu 
adatom |l8| but also for magnetic adatoms in the ma- 
jority spin channel. We also found evidence of an in- 
teraction of the minority spin d-states of the magnetic 
impurities and the surface state continuum, that for Co 
and Ni possibly explains the remarkable enhancement of 
the adatom-induced resonance peak in the minority spin 
channel. The different shape of the adatom-induced reso- 
nance for Cu and Co was indeed observed in experiments 
(l?| and is in qualitative agreement with our calculations. 
Finally, we pointed out the existence of long-range spin- 
polarized oscillations of the surface states around a mag- 
netic impurity being the very origin of a 2D RKKY in- 
teraction. 
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